The design of the experiments (table with parameters) can be found here. But here is a brief description.

PS N_GENES_INITIAL BITING_RATE_MEAN N_LOCI N_ALLELES_INITIAL
01 120 0.1 2 12
02 1200 0.2 2 120
03 12000 0.5 2 1200

I kept daily biting rates fixed to 1.

1. Basic variables for each scenario across regimes

1.1 diversity

In the complete neutrality, the number of circulating strains seems just too much. This is orders of magnitude more than the selection.

Selection

## [1] "PS: 01 | Scenario: S | exp: 001 | run: 1"
## [1] "PS: 02 | Scenario: S | exp: 001 | run: 1"
## [1] "PS: 03 | Scenario: S | exp: 001 | run: 1"

Complete neutrality

## [1] "PS: 01 | Scenario: N | exp: 001 | run: 1"
## [1] "PS: 02 | Scenario: N | exp: 001 | run: 1"
## [1] "PS: 03 | Scenario: N | exp: 001 | run: 1"

Generalized immunity

## [1] "PS: 01 | Scenario: G | exp: 001 | run: 1"
## [1] "PS: 02 | Scenario: G | exp: 001 | run: 1"
## [1] "PS: 03 | Scenario: G | exp: 001 | run: 1"

1.2 EIR

The EIR for selection seems too low? For neutrality it is an order of magnitude more.

Selection

Complete neutrality

Generalized immunity

2. Compare scenarios for a given regime

Notice the order of magnitude difference in the number of repertoires

Low diversity & transmission

## [1] "Scenario: S | exp: 001 | run: 1"
## [1] "Scenario: N | exp: 001 | run: 1"
## [1] "Scenario: G | exp: 001 | run: 1"

Medium diversity & transmission

## [1] "Scenario: S | exp: 001 | run: 1"
## [1] "Scenario: N | exp: 001 | run: 1"
## [1] "Scenario: G | exp: 001 | run: 1"

High diversity & transmission

## [1] "Scenario: S | exp: 001 | run: 1"
## [1] "Scenario: N | exp: 001 | run: 1"
## [1] "Scenario: G | exp: 001 | run: 1"

3. Edge weight values

Selection

Complete Neutrality

For now cannot plot the high diversity because there are just too many edges (computer crashes)

Generalized immunity

4. Modularity

This will be influenced by cutoff. For now I am using the cutoffs I used in the previous version: 0.25 (low), 0.7 (medium), 0.9 (high).

## [1] "04 | S | 001 | 1 | 0.25"
## [1] "05 | S | 001 | 1 | 0.7"
## [1] "06 | S | 001 | 1 | 0.9"
## [1] "04 | N | 001 | 1 | 0.25"
## [1] "05 | N | 001 | 1 | 0.7"
## [1] "06 | N | 001 | 1 | 0.9"
## [1] "04 | G | 001 | 1 | 0.25"
## [1] "05 | G | 001 | 1 | 0.7"
## [1] "06 | G | 001 | 1 | 0.9"

4.1 Example of modules

4.2 Module and strain persistence

Relative persistence

Persistence of repertories is after removing duplicates. Gray is repertoires.

Absolute persistence

Persistence of repertories is after removing duplicates. Gray is repertoires.